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Abstract 

We consider the heat equation in a multidimensional domain with nonlocal hysteresis feedback control 
in a boundary condition. Thermostat is our prototype model. We construct all periodic solutions with 
exactly two switching on the period and study their stability. Coexistence of several periodic solutions with 
different stability properties is proved to be possible. A mechanism of appearance and disappearance of 
periodic solutions is investigated. 
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1 Introduction 

Hysteresis operators arise in mathematical description of various physical processes [HEBUS]. Models with 
hysteresis for ordinary differential equations were considered by many authors (see e.g., [1,2,6, 19, 22, 24l[25]). 
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Partial differential equations with hysteresis have also been actively studied during the last decades (see [4"l[2l>] 
and the references therein). The primary focus has been on the well-posedness of the corresponding problems 
and related issues (existence of solutions, uniqueness, regularity, etc.). However, many questions remain open, 
especially those related to the periodicity and long-time behavior of solutions. 

In this paper, we deal with parabolic problems containing a discontinuous hysteresis operator in the bound- 
ary condition. Such problems describe processes of thermal control arising in chemical reactors and climate 
control systems. The temperature regulation in a domain is performed via heating (or cooling) elements on 
the boundary of the domain. The regime of the heating elements on the boundary is based on the registration 
of thermal sensors inside the domain and obeys a hysteresis law. 

Let v(x, t) denote the temperature at the point x of a bounded domain Q C K n at the moment t. We define 
the mean temperature v(t) by the formula 

v(t) = / m(x)v(x,t) dx, 
Jq 

where m is a given function from the Sobolev space H 1 (Q) (see Condition 1 2. II for another technical assumption 
on m(x)). 

In our prototype model, we assume that the function v(x,t) satisfies the heat equation 

v t (x,t) = Av(x,t) {xeQ,t>0) (1.1) 

and a boundary condition which involves a hysteresis operator % depending on the mean temperature v. 

The hysteresis H(v)(t) is defined as follows (cf. p~9l[26] and the accurate definition and Fig. 12. II in Sec. [2]). 
One fixes two temperature thresholds a and j3 (a < (3). If v(t) < a, then H(v)(t) = 1 (the heating is switched 
on); if v(t) > /3, then H(v)(t) = —1 (the cooling is switched on); if the mean temperature v(t) is between a 
and /3, then H(v)(t) takes the same value as "just before." We say that the hysteresis operator switches when 
it jumps from 1 to —1 or from —1 to 1. The corresponding time moment is called the switching moment. Note 
that the hysteresis phenomenon takes place along with the nonlocal effect caused by averaging of the function 
v{x, t) over Q. 

To be definite, let us assume that one regulates the heat flux through the boundary dQ. Then the boundary 
condition is of the form 

= K(x)H{v)(t) (x e <9Q, t> 0), (1.2) 

ov 

where v is the outward normal to dQ at the point a;, if is a given smooth real- valued function (distribution of 
the heating elements on the boundary). 

A similar mathematical model was originally proposed in 9,10^ Generalizations to various phase-transition 
problems with hysteresis were studied in [4,5,7, 15,20 . Some related issues of optimal control were considered 
in J3J. The most important questions here concern the existence and uniqueness of solutions, the existence of 
periodic solutions, and long-time behavior of solutions. The latter two questions are especially difficult. 

In the case of a one- dimensional domain Q (a finite interval, n = 1), the periodicity was studied in [8[ 
[H1Q3I23] ' Problems with hysteresis on the boundary of a multidimensional domain (n > 2) turn out to be 
much more complicated. Although one can relatively easily prove the existence (and sometimes uniqueness) of 
solutions, the issue of finding periodic solutions is still an open question. The main difficulty here is related to 
the fact that, in general, the solution does not depend on the initial data continuously. The reason is that the 
solution may intersect the "switching" hyperplane {ip = a} or {tp = /?} nontransversally (cf. [2j[25], where the 
same phenomenon occurs for ordinary differential equations) . This leads to discontinuity of the corresponding 
Poincare map. As a result, most methods based on fixed-point theorems do not apply to the Poincare map. 

One possible way to overcome the nontransversality is to consider a continuous model of the hysteresis 
operator. This was done in |12) . where a thermocontrol problem with the Preisach hysteresis operator in 
the boundary condition was considered and the existence of periodic solutions and global attractors were 
established. Note that the periodicity and the long-time behavior of solutions were also studied in [17l[28] in 
the situation where a hysteresis operator enters a parabolic equation itself (see also [27] and the references 
therein) . 

The first results about periodic solutions of thermocontrol problems in multidimensional domains with 
discontinuous hysteresis were obtained in |13j . In [14] . a new approach was proposed. It is based on regarding 
the problem as an infinite-dimensional dynamical system. By using the Fourier method, one can reduce the 
boundary-value problem for the parabolic equation to infinitely many ordinary differential equations, whose 
solutions are coupled with each other via the hysteresis operator. 

In [14] , the existence of a unique periodic solution of the thermocontrol problem is proved for sufficiently 
large (3 — a. This periodic solution possesses certain symmetry, is stable, and is a global attractor. A similar 
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result was established for arbitrary a and ft, but m(x) being close to a constant. The idea was to find an 
invariant region for the corresponding Poincare map and prove that the Poincare map is continuous on that 
region. This turns out to be true for sufficiently large ft — a. However, one can construct examples with small 
ft — a, where an invariant region exists and even is an attracting set, but the Poincare map is not continuous 
on it. 

In the present paper, we will show that the requirement for ft — a to be large enough is essential. We will 
prove that if ft — a is small, then unstable periodic solutions may appear. In particular, they may have a 
saddle structure. To construct those solutions, we will develop a general procedure which yields all periodic 
solutions (with two switchings on the period) in an explicit form. This procedure works even in the presence 
of discontinuity caused by the above nontransversality. In particular, it allows one to find periodic solutions 
on which the Poincare map is discontinuous. 

To study stability of periodic solutions (in particular, to find unstable ones), we propose a method which 
allows one to reduce the original system to an invariant subsystem. The dimension of this subsystem is equal 
to the number of nonvanishing modes in the Fourier decomposition of the m(x). If m{x) has finitely many 
nonvanishing modes, then one can explicitly write down the linearization of the reduced system and find all 
the eigenvalues. They provide complete information about the stability of the periodic solution. 

The invariant subsystem corresponding to the nonvanishing modes of m(x) is called guiding. The remaining 
subsystem is called guided. We prove that the full system (i.e., the original problem) has a periodic solution 
whenever the guiding system has one. Moreover, the periodic solution of the full system is a global attractor 
(is stable, uniformly exponentially stable) whenever the periodic solution of the guiding system possesses 
those properties. We call these results conditional existence of periodic solutions, conditional attractivity, and 
conditional stability, respectively. The above "guiding-guided" decomposition is a result of independent interest. 
It generalizes the results of |13j . where m{x) = const (in our terminology, this corresponds to m{x) which has 
only one nonvanishing mode). 

The paper is organized as follows. In Sec. [2] we define the hysteresis operator, formulate the problem, 
introduce a notion of solution, recall some properties of the solutions, and reduce the problem to an infinite- 
dimensional dynamical system. In the end of Sec. [2j we define the guiding and the guided subsystems and 
introduce the corresponding decomposition of the phase space (the Sobolev space i? 1 (Q)). Most results of this 
section are proved in [14] . 

In Sec. [3l we give a notion of periodic solution with two switchings on the period. By using the Poincare 
maps of the guiding system and the full system, we prove conditional existence of periodic solutions, conditional 
attractivity, and conditional stability. The latter two results are proved under assumption that the periodic 
solution of the guiding system intersects the hyperplanes {(p = a} and {<p = ft} at the switching moments 
transversally. The transversality implies the continuity (and even the Frechet differentiability) of the Poincare 
maps in a neighborhood of the periodic solution. However, we require neither that this neighborhood be 
invariant under the Poincare map, nor that the Poincare map be continuous in a (bigger) invariant neighborhood 
(which exists due to [14]). 

In Sec. IU we show that any periodic solution with two switchings on the period possesses a symmetry in 
the phase space. By using this symmetry, we develop an algorithm which allows us 

1. to construct all periodic solutions (with two switchings on the period) in an explicit form for any given 
a and ft; 

2. to find a sufficient condition under which periodic solutions exist for all sufficiently small ft — a. 

3. to define bifurcation points where periodic solutions may appear or disappear; a role of a bifurcation 
parameter is played either by the period or by the difference ft — a; 

Furthermore, using the results about the guiding-guided decomposition from Sec. [31 we construct examples 
in which periodic solutions are stable or unstable, respectively. In the "unstable" case, we show that they may 
have a saddle structure. 

As a conclusion, we note that the developed method can also be applied to the study of the Dirichlet or 
Robin boundary conditions. Moreover, one can study the problem where the heat flux through the boundary 
(in the case of the Neumann boundary condition) changes continuously. Mathematically, this means that the 
boundary condition (|1.2[) is replaced by 

— = K(x)u(t) (x edQ, t> 0), 
av 

au'(t)+u(t) =H(v)(t) 

with a > (cf. 0[H[T3l[23]). 



3 



2 Setting of the Problem. Reduction to Infinite Dynamical System 
2.1 Setting of the problem 

Let Q C R" (n > 1) be a bounded domain with smooth boundary. Let L2 = Li2(Q)- Denote by H 1 = H 1 (Q) 
the Sobolev space with the norm 



H 1 



(|^)| 2 + |V^)| 2 )^ 



1/2 



Let JJ 1 / 2 = iJ 1 / 2 (9Q) be the space of traces on dQ of the functions from H 1 . 

Consider the sets Q T = Q X (0, T) and r T = dQ x (0, T), T > 0. Fix functions if € ii 1/2 and m G H 1 and 
real numbers a and /3, j3 > a. 

For any function ip(x) or v(x,t) (x € Q, t> 0), the symbol " will refer to the "average" of the function: 



m(x)tp(x) dx, v(t) = / m(x)v(x, t) dx. 
Q JQ 



Let v(x, t) denote the temperature at the point x G Q at the moment t > satisfying the heat equation 

«t(a;, t) = Av(x, t) ((x, t) G Q r ) (2.1) 

with the initial condition 



and the boundary condition 



dv 
dv 



v\ t =o = <p(x) (x G Q) 

= K(x)H(v)(t) ((i,t)er T ). 



(2.2) 



(2.3) 



Here v is the outward normal to Tt at the point (x,t) and is a hysteresis operator, which we now define. 

We denote by BV(0,T) the Banach space of real-valued functions having finite total variation on the 
segment [0, T] and by C r [0,T) the linear space of functions which are continuous on the right in [0,T). We 
introduce the hysteresis operator (cf. [19U26] ) 

% : C[0, T] — > BF(0, T) n C r [0, T) 

by the following rule. For any g G C[0,T], the function /i = H(g) : [0, T] — > {—1, 1} is defined as follows. Let 
X t = {f G (0, t] : g{t') = a or /?}; then 



fe(0) = 



1 if ff(0) < 0, 
-1 if 5 (0)>/3 



and for f G (0, T] 

!/i(0) ifX t = 0, 

1 if ^ and ^(maxXt) = a, 

-1 if Xt 7^ and g(maxX t ) = /3 

(see Fig. 12. A point r such that H(g)(r) ^ %{g){r — 0) is called a switching moment ofH(g). 

n(g) 



- 1 

Figure 2.1: The hysteresis operator % 
We assume throughout that the following condition holds. 
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Condition 2.1. The coefficient K{x) in the boundary condition (j2.3p and the weight function m(x) satisfy 

I K(x)dT>0, j m(x)dx>0. (2.4) 

JdQ JQ 

Remark 2.1. From the physical viewpoint, the function K{x) characterizes the density of the heating (or 
cooling) elements on the boundary and m{x) characterizes the density of thermal sensors in the domain. 
Clearly, inequalities (|2.4j) hold in the physically relevant case K(x) > for a.e. x £ dQ, K(x) ^ 0, and 
m(x) > for a.e. x £ Q, m(x) ^ 0. 

2.2 Functional spaces and the solvability of the problem 

For any Banach space B, denote by C([o,6]; B) (a < b) the space of B-valued functions continuous on the 
segment [a, b] with the norm 

\\u\\ C ([a,b];B) = max \\u(t)\\ B 
and by ^((ctj B) the space of L2-integrable B- valued functions with the norm 

' ,6 \ 1/2 



IMU a ((a,6);B) = \\u(t)\\ B dt J 

We introduce the anisotropic Sobolev space H 2,1 (Q x (a, b)) = {v £ Lzifa, b); H 2 ) : v t £ L2((a, b); L2)} with 
the norm 

\ 1/2 



IMI*p.i(Qx(a,&)) = yj \H-,t)\\ H 2dt + J IM-^)IL 2 dt J 

Taking into account the results of the interpolation theory (see, e.g., [2TJ Chap. 1, Sees. 1-3, 9], we make 
the following remarks. 

Remark 2.2. The continuous embedding H 2,1 (Q x (a, b)) C C([a,b], H 1 ) takes place. Furthermore, for any 
v £ _ff 2,1 (Q x (a, b)) and r £ [a,b], the trace v\ t=T £ H 1 is well defined and is a bounded operator from 
H 2 ^(Q x (a, b)) to H 1 . 

Remark 2.3. Consider two functions v% £ H 2,1 (Q x (a, b)) and v-i £ iJ 2,1 (<5 x (b, c)), where a < 6 < c. 
Let v(-,t) = «i(-,i) for £ G (a, 6) and v(-,t) = V2(-,t) for t £ (b,c). Then w e H 2 ' X {Q x (a, c)) if and only if 

«l|t=6 = W2|i=6- 

Definition 2.1. A function v(x, t) is called a solution of problem (|2.ip ~ (|2.3p in Qt with the initial data ip £ H 1 
if v £ H 2,1 (Qt) and v satisfies Eq. (|2.1I) a.e. in Qt and conditions (|2.2p . (|2 . 3[) in the sense of traces. 

Definition 2.2. We say that v(x, t) (t > 0) is a solution of problem (|2.1l) - (|2.3p m Qoo if it is a solution in Qt 
for all T > 0. 

The following result about the solvability of problem (|2.ip - (|2.3l) is proved in JT4] Theorem 2.2]. 

Theorem 2.1. Let <p £ H 1 and \\cp\\ H i < R (R > is arbitrary). Then there exists a unique solution v of 
problem (12. 1[) — (|2.3[) in Qoq and the following holds for any T > 0. 

1. One has 

|K-,t)|| ff i < c \\v\\ h ,. HQt) < cidl^lffi + \\K\\ Hl/2 ), (2.5) 
where cq = Cq(T) > and ci = Ci(T) > do not depend on ip and R; 

2. The interval (0,T] contains no more than finitely many switching moments t\ < ti < . . . < tj ofH(v). 
Moreover, 

U-U-l <t* + 2(/? ~ a) ' i = l,2,..., (2.6) 
m K 

where t* depends on m and R but does not depend on ip, T, a, f3\ 

, , ^ * ■ J 1,2,..., J if<p<aorip>P, 
[2, 3,..., J ifa<ip<(3 7 

where 

T *= const (2.8) 
wi£/i const > depending on R rather than on m,ip,T,a, /3; in ()2.6|) and (|2.8p . to = 0. 
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2.3 Reduction to infinite-dimensional dynamical system 



Due to Theorem 12.11 the study of the solutions of problem (12.1I) ~ (|2.3I) with hysteresis can be reduced to the 
study of the solutions of parabolic problems without hysteresis by considering the time intervals between the 
switching moments ij. 

Thus, if H{v){t) = 1, then problem p7T ]) -([2T5 ]) takes the form 

v t {x,t) = Av(x,t) ((x,t)eQ T ), (2.9) 

v(x,0) = <p(x) (xeQ), (2.10) 

dv 



dv 



r 7 



K{x) ((i,t)6r T ). (2.11) 



HH{v)(t) = -1, one should replace K(x) by -K(x) in (|2~TT1) . 

Definition 2.3. A function v(x,t) is called a solution of problem (|2.9j) - (|2.11j) in Qt if V £ H 2,1 {Qt) and v 
satisfies Eq. (|2.9p a.e. in Qt and conditions (|2.10j) . (|2.11|) in the sense of traces. 

It is well known that there is a unique solution v € H 2,1 {Qt) of problem (|2.9[) - (12.11l) . 

Now we give a convenient representation of solutions of problem (|2.9p - (|2.11j) in terms of the Fourier series 
with respect to the eigenfunctions of the Laplacian. 

Let {Xj}JL and {ej(x)}°°^ denote the sequence of eigenvalues and the corresponding system of real- valued 
eigenfunctions (infinitely differentiable in Q) of the spectral problem 

dc 

- Aedx) = Aj-e,-(x) (x £ Q), =0. (2.12) 

dv dQ 

It is well known that = Ao < Ai < A2 < . . . < Xj < . . . , eo(x) = (mesQ) -1 / 2 > 0, and the system of 
eigenfunctions {ej}°^ can be chosen to form an ortho normal basis for L^- Then, the functions ej/y/Xj + 1 
form an orthonormal basis for H 1 . 

Remark 2.4. In what follows, we will use the well-known asymptotics for the eigenvalues Xj = Lj 2 / n + o(j 2//n ) 
as j —> +00 (L > and n is the dimension of Q). 

Any function ip £ L2 can be expanded into the Fourier series with respect to ej(x), which converges in L2: 

00 00 
iP(x) =J2^j(x), ML = E l^'l 2 ' ( 2 - 13 ) 

where ipj = f Q ip(x)ej(x) dx. If ip £ H 1 , then the first series in (I2.13[) converges to tp in H 1 and 



3=0 

Denote 

m(x)ej(x) dx, Kj = / K{x)e 3 {x)dx (j = 0, 1,2, . . .). (2.15) 

Note that mo,Ko > due to Condition 12.11 We also note that Kj are not the Fourier coefficients of K(x). 
However, the following is proved in |14j : 

£(W + kW Phv „ ( , 16) 

where c > does not depend on if. 

Remark 2.5. Using (|2.14j) and (|2.16|) . we obtain the estimate 

00 

Eh^l <c\\ m \\ m \\K\\ H1/2 , 
i=i 

which will often be used later on. 
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The numbers rrij and Kj play an essential role when one describes the thermocontrol problem in terms of 
an infinite-dimensional dynamical system. The following result is true (see |14l Lemma 2.2]). 

Theorem 2.2. Let ip 6 H . Then the following assertions hold. 

1. The solution v of problem (|2.9p ~ (|2.1ip can be represented as the series 

(2.17) 



v(x,t) = ~S~^Vj(t)ej(x), t > 0, 
j=o 

where Vj(t) = §Qv(x,t)ej(x) dx andvjit) satisfy the Cauchy problem 

Vj(t) = -Xjvjit) + Kj, Vj (0) = Vj C = d/dt, j = 0, 1, 2 . . .)• 
The series in (|2.17p converges in H 1 for all t > 0. 
2. The mean temperature v(t) is represented by the absolutely convergent series 

oo 

which is continuously differentiate for t > 0. 
Remark 2.6. In what follows, we will also use the explicit formulas for the solutions of Eqs. (|2.18[) 



(2.18) 



(2.19) 



K 



K, 



«0(i) = ip + K Q t, Vj (t) = i^p 3 - j^j e-^ 1 + ^L, i = 1,2,.. 
Formally, relations (|2.18l) can be obtained by multiplying (|2.9[) by ej(x), integrating by parts over Q, and 

oo 

substituting v(x, t) = ^2 Vj(t)ej(x). The rigorous proof is given in [14] . 

3=0 

A geometrical interpretation of the dynamics of vo(t),vi(i), ... is as follows. We choose the orthonormal 
basis in L2 (which is orthogonal in H 1 ) consisting of the eigenfunctions eo, e\, €2, • • ■ ■ Then, in the coordinate 
form, we have 

e = (1,0,0,0,...), e x = (0,1,0,0,...), e 2 = (0,0,1,0,...), ... 

and (cf. IpTf]) ) 

<p = (tpo,<p 1} <p2,...), v(-,t) = (v (t),v 1 (t),v 2 (t),...). 

Consider the plane going through the origin and spanned by the vector eo = (1, 0, 0, ... ) and the vector 
m = (mo, mi, m2, ■ ■ ■ ) (if they are parallel, i.e., mi = TO2 = • •• =0, then we consider an arbitrary plane 

containing eo). We note that the angle between the vectors m and eo is acute (their scalar product is equal to 

00 

mo > 0). Clearly, the orthogonal projection of the hyperspace <p — ^2 m j L Pj = a ( or P) on this plane is a line 

3=0 

(see Fig. 

Due to (|2.18p . vo(t) "goes" from the left to the right with the constant speed Kq > 0, while Vj(t) exponen- 
tially converge to Kj/Xj (see Fig. I2.3|) . 



v(.,t) 




tp = a 



Figure 2.2: The plane spanned by eo = (1,0,0,...) and m = (mo, mi, ma, . . . ). 

Due to Theorem 12.21 the original problem (|2.ip ~ (|2.3p can be written as follows: 

v Q (t)=H(v)(t)K Q , v (0) = <p , 

Vj (t) = -X jVj (t) + U{v)(t)K v Vj (0) = tpj (j = 1, 2 . . . ). 

Equations (|2.20p define an infinite-dimensional dynamical system for the functions Vj (t) . These functions are 
"coupled" via formula (|2.19l) for the mean temperature, which is the argument of the hysteresis operator H. 



(2.20) 
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Figure 2.3: The plane spanned by e, and e 3 -, i ^ j, i, j > 1. 
2.4 Invariant subsystem and "guiding-guided" decomposition 

In this subsection, we show that if some coefficients nij vanish, then the system (|2.20[) has an invariant 
subsystem. 

We introduce the sets of indices 

J = {j e N : mj ^ 0}, J = { j e N : = 0}. 

Clearly {0}UJU Jo = {0,1,2,...}. 

Note that, for any solution v(x,t) of problem (12.1l) - (|2.3p . we have (cf. (J2TT9J) ) 

o(t)= E ™,^(i), *>o. 

je{o}uj 

Therefore, the dynamics of Vj(t), j € {0} U J, does not depend on the functions Vj(t), j € Jo, and is described 
by the invariant dynamical system 

v {t)=H{v){t)Ko, v (0) = <p , 

v j (t) = -\ j v j (t)+H(W)Kj> «i(o) = *>j (ieJ). 

The dynamics of Vj(t), j € Jo, is described by the system 

w i (t) = -A i i; i (t)+W(fi)(i)lf J -, «i(0) = vj (ieJo), (2.22) 
where the hysteresis operator % depends only on the functions Vjit) from the system (|2.21|) . 
Definition 2.4. We say that the system (|2.21|) is guiding, while the system (|2.22|) is guided (by (|2.2ip V 

In what follows, we will use the following notation. For any number ipo and (possibly, infinite-dimensional) 
vectors {fj}jel an d {<Pj}je3 (<Pj we denote 

<P = {<Pj}jel> H> = Wj}je{o}uh Vo = {^iljeJo- 

Thus, e.g., v(i) and Vo(t) will represent the solutions of the guiding system (|2.21[) and the guided system (|2.22[) . 
respectively. 

The above decomposition of the system (|2.20|) implies the corresponding decomposition of the phase 
space H x : 

H 1 = R x V x Vo = V x V Q , (2.23) 
where the norms in V, V, and Vo are given by 

I « 1/2 / v 1/2 / N 1/2 

\\v\\v= (E( 1+ ^)i^i 2 ) . n*% = ( E (i+^ii 2 ) > iivoiivb = (E( 1+A i)i^i 2 

\ieJ / \ie{o}uj / \jeJo 

(2.24) 

Further we show that Definition 12.41 is quite natural. In particular, we prove that if z(t) is a periodic 
solution of the guiding system (|2.21[) . then there exists a periodic solution of the full system (|2.20[) of the 
form (z(i), zo(t)). Moreover, the latter is stable if and only if z(t) is a stable periodic solution of the guiding 
system (pHiTj) . 

As an application of this result, assuming that the set J is finite, we will construct a periodic solution z(x, t) 
with small period such that z(t) is an unstable periodic solution of the guiding system (|2.21l) . Clearly, z(x,t) 
will be unstable in this case, too. 
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3 Periodic Solutions 



3.1 Conditional existence of periodic solution 

We begin with a definition of periodic solutions (with two switchings on the period) of problem (j2.ip . (|2.3[) . 
Recall that the symbol " refers to the "average" of the function (see Sec. 12. ip . 

Definition 3.1. A function z(x, t) is called an (s, a)-periodic solution (with period T = s + a) of problem (|2.1I) . 
(|2.3|) if there is a function ip g H 1 such that the following holds: 

1. ip = a, 

2. z(x,t) is a solution of problem (I2.1j) - (I2.3[) (in Qoo) with the initial data ip, 

3. there are exactly two switching moments s and T of H(z) on the interval (0, T] (such that z(s) = f3 and 
z(T) = a), 

4. z(x,T) = z(x,0) (= ip(x)). 

Definition 3.2. If z(x, t) is an (s, cr)-periodic solution of problem (|2.ip . (|2.3p and T = s + <r, then the sets 

T = {z(; t), t e [0, T]}, f = {z(t) : t € [0, T]}, T = {*„(*) : * € [0, T]}. 
are called the trajectories of z(x,t), z(t), and Zo(i), respectively. 

We also consider two parts of the trajectory corresponding to the hysteresis value W(z) = 1 and — 1: 

Y X = {z(;t), t e [0, *]}, T 2 = {z(; t), * G [ S , T]}. 

Similarly, one introduces the sets and Toj, j = 1,2. 

Remark 3.1. It follows from the definition of the hysteresis operator H and from Definition 13. II that if z(x, t) 
is an (s, <r)-periodic solution with period T = s + a of problem (|2.ip . (|2.3p . then 

H(2)(t) = l, te[0,s); K(i)(t) = -1, te[s,T). 

Further in this section, we establish the connection between periodic solutions of the guiding system (|2.21l) 
and those of the full system (|2.20[) . The definitions of (s, er)-periodic solutions for the guiding system (|2.21l) 
and for the full system (|2.20l) are analogous to Definition 13.11 

The following theorem generalizes Theorem 4.4 in [T3], where mo 7^ and mi = m 2 = ■ ■ • = (i.e., J = 
and Jo = N). 

Theorem 3.1. Let z(t) be an (s, a) -periodic solution of the guiding system (|2.21j) . Then there exists a unique 
function Zo(t) such that (z(t),Zo(t)) is an (s, a) -periodic solution of the full system (|2.20p (which generates an 
(s, a) -periodic solution z(x,t) of problem (|2.ip . (I2.3P ). 

Proof. We recall that the spaces V and Vq form the decomposition of H 1 (cf. (|2.23p ). 

We introduce a nonlinear operator Mt : Vq — > Vq as follows. For any <p e Vq, we consider the element 
(z(0), y>o) G -ff 1 - By Theorem 12.11 and the invariance of the guiding system (|2.2ip . there is a unique solution 
of the full system (|2.20p . which is of the form (z(£),vo(t)) € H 1 . Clearly, Vo(i) is a solution of the guided 
system (|2~2"2"|) . We set 

Mt( Vq ) = MT), T=s + a. 

We claim that Mf is a contraction map. Indeed, let <pj, <Po G ^0 and let (i) and Vq (i) be the corresponding 
solutions of the guided system (|2 . 22|) . Since the mean temperature is defined via z(t) and does not depend on 
Vq(£) and Vq(£), it follows that the difference Wo(i) = vj(i) — Vg(i) satisfies the equations 

wj(t) = -\jWj(t), Wj(0) = <p] - (ieJo). 

Therefore, 

||M T ( Vo ) - M T (^)|| 2 Vo = ||wo(T)|| 2 yo = 2(1 + A J )e- aA ^|^ - ^ 2 | 2 < e- 2 ^|| Vo - ^ll^ , 

where >f = minA, > 0. 

jeh 

Thus, has a unique fixed point ip € Vq, which yields the desired (s, cr)-periodic solution (z(t) , Zo(t)) 
of the full system (j2~20l) . □ 

Further, we will study the connection between the stability and attractivity of solutions of the guiding 
system and the guided and full systems. To do so, we need to define the Poincare maps of the respective 
systems. 
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3.2 The Poincare maps 

In this subsection, we introduce the Poincare maps for the full system (|2.20|) and for the guiding system (|2.22[) . 
It is proved in [14] that the stability of a periodic solution of the full system follows from the stability of the 
corresponding fixed point of the Poincare map. Therefore, we will concentrate on the properties of the Poincare 
map. 

We consider nonlinear operators (see Fig. 13. ip 

P a : {ip e H 1 : <p < (3} ^ {<p e H 1 : <p = f3}, 
Pp : {cp G H 1 : <p > a} -> {cp G H 1 : Cp = a} 

defined as follows. 

Let cp G H , cp < (3, and let v(x,t) be the corresponding solution of problem (|2.1|) - (|2.3|) in (Qoo)- Due to 
Theorem 12.11 there exists the first switching moment t\ such that v(t\) — j3 and there are no other switchings 
on the interval (0, t\). In other words, the function v a (x, t) := v(x, t) is a solution of the initial boundary-value 
problem on the interval (0,ii): 

v?(x,t) = Av a (x,t) ((x,t)eQ tl ), (3.1) 
v a (x,0) = cp{x) (xgQ), (3.2) 

^=K(x) ((Mler,). (3.3) 

We set P a {cp) = v a (-,h). 

The operator is defined in a similar way. Let cp G H 1 and cp > a. As before, there is a moment t 2 > 
and a function v@(x,t) such that v^{x,t) is a solution of the problem 



vP{x,t) = AvP(x,t) ((x,t)eQ T2 ), 

v l3 (x,0) = cp(x) (xeQ), 
dv 13 

— = -K{x) i,fer„ 

ov 

vP{ji) > a for t < r 2 , and -W 3 (r 2 ) = a. We set Pp(cp) = w /3 (-,t 2 ). 

We introduce the Poincare map for problem (|2.1[) - (|2.3[) . or, equivalently, for the full system (12.20[) 

P : {cp G H 1 : <p < /3} -> {ip G H 1 : cp = a}, 
P = PpP a . 



(3.4) 
(3.5) 

(3.6) 



ei,e 2 , 



t 


k ^ — ' 
















. X 

cp = a \ m 


= 





{Ki/^KjlXj) 




{-Ki/^-Kj/Xj) 

Figure 3.1: The operators P Q and P = P/jPa on the planes (eo, m) and (e^, ej), i =/= j 



We also introduce the operator (functional) ti : {ip G H 1 : cp < (3} — > M. given by 

ti(cp) = the first switching moment of T-L{v) for system (|2.ip - (j2.3[) . 
We will use the following result (see Remark 4.3 in [14]). 
Lemma 3.1. Let z(x,t) be aperiodic solution of problem (|2.1[) . (|2.3I) . // 

dz 

— ^ at the switching moments, 

then the operators P a (cp), P(cp), and ti(cp) are continuously differentiable in a neighborhood ofT\ n {(p G H 1 
cp < 0}. The operator P p{cp) is continuously differentiable in a neighborhood o/T 2 D {cp E H 1 : cp > a}. 
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We denote by E : H 1 — » V the orthogonal projector from H 1 onto V. 

Similarly to the operators P a , P^, P, and ti, we introduce the operators n Q , LT^, II, and ti, respectively, 
corresponding to the invariant guiding system ()2.21|) and defined on the elements from V. Due to the invariance 
of we have 

U a (<p) = EP a (<p,<p ), n fj (<p)=EP l3 ((p,<p ), fi(£)=EP(£, Vo ), ti(<p) = ti(ip,<p ) 

V£> e v, vv e v . 

We say that II is the guiding Poincare map. 

The following theorem shows that the stability of (exponential convergence to) a fixed point of the guiding 
Poincare map II implies the stability of (exponential convergence to) the fixed point of the Poincare map P of 
the full system. 

First, we introduce some notation. Let z(x,t) be an (s, cr)-periodic solution with period T = s + a of 
problem (|2.ip . (|2.3|) and (z(t),Zo(t)) the corresponding periodic solution of the full system (|2.20p . We denote 

t/> = z(0), Vo=zq(0). 

Let u(a;, t) be another solution of problem (|2.1I) - (|2.3I) such that -0(0) = a, and let (v(i),vo(t)) be the 
corresponding solution of the full system (|2.20[) . The initial data will be denoted by 

¥> = v(0), <p = v (0) 

and the consecutive switching moments by t\,t%,.... We also set £ = 0. 
Theorem 3.2. Suppose that 

dz 

— 7^ at the switching moments. 

1. For any 8$ > 0, there exists S > such that ifv(ti) remain in the 5 -neighborhood of xj) for even i and in 
the 8 -neighborhood o/z(s) for odd i [i = 0, 1, 2, . . . ), then, for all ip Q in the 5q -neighborhood of ip , Vo(£j) 
remain in the 8^-neighborhood of ij) for even i and in the 5q -neighborhood of zq(s) for odd i; 

2. Let 

\\^{U)-^\\v + M^i)-^s)\\^<kt, i = 0,2,4..., (3.7) 

for some < q < 1 and k > which do not depend on i. Then, for any neighborhood Vq of tj> and for 
all (p € Vo , 

\\v (ti)-il> \\ Vo + \\v (ti+i)-z{s)\\v <koqi, i = 0,2,4..., (3.8) 

where < qo — qo(q) < 1 a^rf fco = fco(fc: Vo,q) > do not depend on (p and ip in the corresponding 
neighborhoods. 

In the proof of this theorem, we will use the following technical lemma. 
Lemma 3.2. Let a sequence bo, bi, &4, • • • of nonnegative numbers satisfies the inequalities 

bi + 2<tbi + kv\ i = 0,2,4,..., 

where k > and < C? ^ < 1 do not depend on i. Then there are numbers < q = q((,,v) < 1 and 
c = c(C, v, q) > which do not depend on i such that 

b i <(b + c)q i , i = 0,2,4,.... 

Proof. Let 7 = max(C, ^ 2 ). Clearly, < 7 < 1. Then 

fo 4+2 < ih + kf' 2 , i = 0,2,4,.... 

Consider the sequence Cj = bi^ 1 / 2 . It satisfies 

Ci+2 < Q + fc7~ 1 < co + fcii, 

where k\ > does not depend on Ci and i, which yields the desired estimate of 6j. □ 
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Proof of Theorem \3.S\ 1. Let us prove assertion 1. 

la. By assumption, s = ti(i/>) is the first switching moment of H(z). Denote by t = t\ = ti(<p) the 
first switching moment of H(v). By Lemma [3.11 ti(<p) is continuously differentiable in a sufficiently small 
(^-neighborhood of if). Hence, 

\r - s\ < ki\\<p - 4>\\y < hS, 



(3.9) 



where ki > depends on 5\ but does not depend on 8 < S±. 

lb. Using Remark 12.61 and the fact that, for any e > 0, there is d e > such that 



{a + bf < (1 + e)a 2 + d e b 2 



we obtain 



|v (r) - M')\\v = E( X + X i)M T ) z ^ 



E( 1 + A ^ 



K 
% 



3 » e" A ^ 



A' 
A^ 



(3.10) 



< (1 + e) + A,)|^- - e" 2A ^ + 4 E(l + Aj) 
jel jeJo 



Now we fix e > such that 



C= (l + 2e)e- 2 " s < 1, 



where x = min Xj > 0. Further, taking into account (|3.9|) . we choose 5 > so small that 

(l + e)e- 2 ^ T < (l + 2e) e - 2 ^ s . 
Combining (j3~T0| . ([3~TH) . and (|3~T2|) yields 



||v (r) - as (*)||y < CIIVo - ^ollvo + ^ E (! + A J 
Using ([2~T6]) . (j2~24l) . and estimate (JUSJ), we deduce from (|37T3|) 



1>i - — 



A, 



,-A; 



(3.11) 



(3.12) 



(3.13) 



||v (r) - z (s)\\ 2 Vo <(S 2 +k 2 (S), 

where ^((5) > and fe(^) — > as (5 — > 0. In particular, this implies that Vo(ti) = Vo(r) belongs to the 
Jo-neighborhood of zo(s), provided that (5 = <5(<5o) is sufficiently small. 

In the same way, one can now show that Vo(ta) belongs to the (^-neighborhood ofip = Zq(T). By induction, 
we obtain assertion 1. 

2. Now we prove that 

||v (ii)-Volk < *0«o s * = 0,2,4.... (3.14) 

The rest part of estimate (|3.8p can be proved analogously. 

2a. First, we assume that v(0) is in a sufficiently small ^-neighborhood of if). Then, similarly to p,13p . we 
have for even i 



\v (t i+1 ) - Z (s)\\ 2 Vo < C\\MU) - Mv + ^ 53( X + A J 

jela 



^ J A, 



- An S I 



where r»+i = £i+i — U = ti(v(i,)). Using the differentiability of ti and estimate (13.71) . we have 

\n+i - s\ < fei||v(ti) - -0H ^ < fc 3 <f 
Due to (|3.16p . we can assume that rj+x > s/2. Then, taking into account (|3.16p . we have 



(3.15) 



(3.16) 



-AtfT*. 



"An S 



< Aje" A ^/ 2 |T 4+1 -s\ < k. 



Combining this inequality with (|3T5]l . (|2l23|) . and (l2~T6|) yields 

||v (i t+1 ) - z ( S )|| 2 , < Cl|v (*i) - ^ollv& + 
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Making one more step and using the last inequality, we obtain 

||v (ti +2 ) - Volly„ < Cl|v (t t+ i) - z ( S )|| 2 v . + k 5 q 2 * < aCWMU) - ^ollvo + fc 59 2i ) + k 5 q 2i 
<C 2 \\Mti)-^\\ 2 Va +k 6 q 2 \ 

Due to Lemma \'3. 21 the latter inequality implies (|3.14[) . 

2b. Now we take an arbitrary v(0) in the ^-neighborhood of xj). Due to (|3.7[) . there exists an even number / 
(which does not depend on v(0)) such that v(tj) is in the ^-neighborhood of %]). Then the inequality in (|3.14|) 
holds for i = I, I + 2, 1 + 4, . . . due to part 2a of the proof. 

Theorem 12.11 implies the existence of 6 > (which depends on k and Vo but does not depend on v(0) and 
vq(0)) such that ti < 9. Furthermore, Theorem 12.11 implies that 

max ||v (t)||v < *7(M.Vo) = *b(*,Vo). 

tG [0.6] 

Hence, the inequality in p,14|) holds for i = 0, 2, 4, □ 

Remark 3.2. It follows from the proof of Lemma 13.21 that the convergence rate q is greater than 7 1 / 2 = 
max(C 1 / 2 , v) but can be chosen arbitrarily close to this number. 

Therefore, the convergence rate qo in estimate Q3.8P is greater than max (e _> ", e~'*( T ~ s \ qj but can be 
chosen arbitrarily close to this number. 

3.3 Conditional attraction and stability of periodic solution 

Let z(x,t) and v(x,t) be the same as above, but now we do not assume that w(0) is necessarily equal to a. 

The following theorem shows that the convergence to the periodic orbit in the guiding system implies the 
convergence to the corresponding periodic orbit in the full system. Thus, we call the phenomenon in that 
theorem the conditional attraction. 

For the trajectories, we will use the notation given in Sec. 13.11 

Theorem 3.3. Suppose that 



Let 



dz 

— i= at the suntchmq moments, 
dt r y 

fdist(v(t),fO<fc^ if H(v)(t) = l, 
\dist(v(t),f 2 )<fc$ 4 if H{v)(t) = -1, 

for some < q < 1 and k > 0. Then, for any bounded set Vo in Vq, there exist < q — q(q) < 1 and 
k = k(k, Vo, q) > such that, for all ip Q G Vo and t > 0, 

fdist(»(-,i),ri) < kq* if H(v)(t) = l, 
\dist(«(-,t),r 2 )< kq* if H(v)(t) = -1, 

Proof. 1. Suppose we have shown that 

\\v(-,U)-z{-,0)\\ H i + \\v(;t i+1 )-z(;s)\\ H i < kq\ i = 0, 2, 4 . . . , (3.18) 

where < q — q(q) < 1 and k = k(k,Vo,q) > 0. Then, using Lemma 13.11 and arguing as in the proof of 
Theorem 4.3 in |14) . we complete the proof. 
So, let us prove estimate (|3. 181) . 

2. Consider the intersection of the closure of Ti with the set {ip £ V : <p — /?}. This intersection consists of 
the single point z(s), where s is the switching moment of the periodic solution. 

dz 

Since — ^ 0, it follows from the implicit function theorem that there exist a number L > and a 

dt t=s 

sufficiently small number do > such that if 

\/3-z(r)\<d 

for some d < do and r € [0, s], then 

It - s\ < Ld. 
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3. Consider i = 1, 3, 5, ... . By assumption, there exists r € [0, s] such that 

\\v(U)-Z(t)\\v <~kq u . (3.19) 

This inequality together with the Cauchy-Bunyakovskii inequality implies that there exists a constant C\ > 
such that 

|/3 - l(r)| = |0(t,) - z(r)| < CM U - (3.20) 

3a. First, we assume that ti > 9, where 9 > is so large that C\kq s < do (do is the number from part 2 of 
the proof). Then, due to part 2 of the proof, we have 

\t - s\ < LC\kq u . (3.21) 

It was proved in [Ml Lemma 4.6] that the periodic solution z(t) is uniformly Lipschitz-continuous on [0,T], 
which (together with Q3.2ip ) implies that 



z(t) - z(s)\\y < C 2 kq u (3.22) 



for some C 2 > 0. 

Estimates (|3TT^|) and yield 



\\v(U) - z(s)\\y < k(l + C 2 )~q u for U > 9. 
Taking into account (|2.7I) . we have 

\\v(U) - z(s)\\ v < hq\ fovU>e, (3.23) 

where k\ > k and < q\ < 1. 
3b. For ti < 9, we have 

\\v(ti) - z(s)\\ v < dist(v(ti), f x ) < k. (3.24) 
Combining (f3723]) and (1334)) yields 

||v(ti) - z(s)||^ < fe 2 Qi for alltj. 

Applying similar arguments to v(ti+i) and z(0) and using Theorem 13.21 we obtain (|3.18[) . □ 

Now we discuss the phenomenon of conditional stability. When studying the stability of periodic solutions, 
one considers its small neighborhood. When doing so, one has to take into account the initial state of the 
hysteresis operator. 

Definition 3.3. An (s, <r)-periodic solution z(x,t) of problem (|2.1j) . (|2.3j) is stable if, for any neighborhoods 
U\ of Ti and U2 of T 2 in H 1 , there exist neighborhoods Vi of Ti and V2 of T 2 in H 1 such that if 

tp G Vi, if < j3 or cp e V 2 , (p > /3, 

then the solution v(x,t) of problem (|2.ip - (|2.3p in Qrx, with the initial data cp satisfies for all t > 0: 

fwe^i if W(«)(*) = l, 
|weW 2 if H(v)(t) = -l. 

An (s, er)-periodic solution is unstable if it is not stable. 

Definition 3.4. An (s, cr)-periodic solution z(x,t) of problem (|2.ip . (|2.3|) is uniformly exponentially stable if 
it is stable and there exist neighborhoods VVi of Ti and W2 of T 2 in i/ 1 and numbers < q < 1 and fc > 
such that if 

ipeyVi,tp</3 or f£W 2l (p>!3 1 
then the solution v(x,t) of problem (|2.1I) - (|2.3I) in Qao with the initial data if satisfies 

fdist(t>(-,t),ri) < kq l if H(v)(t) = l, 
\dist(v(;t),T 2 ) < kq l if H{v)(t)=-1 

for all t > uniformly with respect to if. 
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Let z(t) be an (s, cr)-periodic solution of the guiding system (|2.2ip . Then, by Theorem 13. 1[ there exists a 
unique function Zg(t) such that (z(t),Zo(t)) is an (s, tr)-periodic solution of the full system (|2.20j) . We denote 
by z(x,t) the corresponding (s, <r)-periodic solution of problem (12.11) . (|2.3[) . 

Theorem 3.4. Suppose that 

dz 

— f= at the switching moments. 
Then the following assertions are equivalent. 

1. The periodic solution z(x,t) of problem (|2.1|) . ()2.3|) is stable (uniformly exponentially stable). 

2. The periodic solution z(t) of the guiding system (|2.21[) zs stable (uniformly exponentially stable). 

3. The element z(0) is a stable (uniformly exponentially stable) fixed point of the Poincare map II. 

Proof. Implication 1 2 is obvious. 

Implication 2 => 3 is proved similarly to the proof of Theorem 13.31 

To prove implication 3 4 1, one should use Lemma 13.11 and Theorem 13.21 and argue as in the proof of 
Lemma 4.7 and Theorem 4.4 in [14] . □ 



4 Symmetric Periodic Solutions 
4.1 Preliminary considerations 

It was noted in [14] that any (s, cr)-periodic solution possesses a certain symmetry, provided that it is unique. 
In fact a much stronger result holds, namely, we show that any (s, (r)-periodic solution possesses symmetry. 

We underline that the results in the previous sections did not depend on the symmetry of periodic solutions, 
but the results of this section do. In particular, by exploiting the symmetry, we give an algorithm for finding 
all periodic solutions with two switchings on the period. Using their explicit form, we will study their stability. 

Definition 4.1. An (s, er)-periodic solution z(x, t) of problem (|2.1[) . (|2.3j) is called symmetric if zj(0) = —Zj(s), 
.7 = 1,2,.... 

Lemma 4.1. Let z(x,t) be an (s, a) -periodic solution of problem (12.11) . (|2.3p . Then s — a and z(x,t) is 
symmetric. 

Proof. Let ip(x) — z(x, 0) = z(x, s + a) and £(x) = z(x, s). 
By Remark 12.61 

£ =ipo + K s, (4.1) 

i, (<% ^ • ^- 3 1- (4-2) 

Applying Remark l2.6l (with Kj replaced by —Kj), we conclude that 

V>o = Co - K a, (4.3) 

o • x)' A 7 x- 1 '• (4 - 4) 

Equalities (|4.ip and (|4. 3[) imply that s = a. Summing up (|4.2j) and (14. 4p and taking into account that s = a, 
we see that 

^• + ^- = (% + 0) e ^ s ; i>i- 

Hence, ^ = — and z(a;,i) is symmetric. □ 

Remark 4.1. Lemma l4Tl shows that the period (and the second switching time) of any (s, cr)-periodic solution 
is uniquely determined by the first switching time and vice versa. Therefore, we will say "I '-periodic solution" 
or just "periodic solution" instead of saying "symmetric (s, s)-periodic solution with period T = 2s" . 

In |14j . it was shown that there is a number 8i > such that if (3 — a > Si, then there exists a periodic 
solution of problem (|2.ip . (|2.3I) . Furthermore, there is a number 62 > Si such that if f3 — a > 62, then there 
exists a unique periodic solution of problem (|2.1j) . ()2 .3|) : moreover, it is stable, and is a global attractor. Both 
numbers £1 and S2 depend on Q, m, and K. 

In this section, we will formulate a sufficient condition which may hold for arbitrarily small j3 — a and still 
provides the existence of (symmetric) periodic solutions. We will show that these solutions may be both stable 
and unstable. 
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Lemma 4.2. Let z{x,t) be a solution of problem (I2.1[) - (|2.3[) with the initial data ip, ip = a, and let s > be 
the first switching moment ofH(z). If 

^j( s ) = -^j, J = 1, 2, ... , 
then z{x,t) is a {symmetric) 2s-periodic solution of problem (12.1[) . (|2.3[) . 

Proof. 1. First, we show that there are no switchings for t G (s, 2s) and that the second switching occurs 
exactly for t = 2s. To do so, we have to show that z(t) > a, or, equivalently, z(s) — z(t) < (3 — a for t E (s, 2s). 
Using Remark 12.61 (with Kj replaced by —Kj) and the assumption that Zj(s) = —ipj, we have for t € (s,2s) 



Zj(t) = Zj (s) + 



El 

A, 



A, 



-fa 



K, 



-Aj(t-s) 



Kj 



J = 1,2,.. 



(4.5) 



«b(t) = z (s) - K [t - s) = ip + 2K s - K t. 
Therefore, taking into account (|2.19|) . we have 



i(s) - z(t) = m K {t - s) + ^2 m j ( $3 ~ y 1 ) ( 

j=o ^ 3 ' 



m^K^aQ + irij ( ip 



3=0 



El 

A, 



e A ^*- s ) -l) 
(e-^ e - 1) = - i>, 



(4.6) 



where 9 = t — s € (0, s). But 5(0) — -0 < (3 — a for 8 £ (0, s) and z(s) — t/j = /3 — a (because s is the first 
switching moment by assumption). 

2. Now we show that z(x, 2s) = ip(x). Indeed, using (|4.5I) and the assumption that Zj(s) = —ipj, we obtain 



z 3 (2s) = - 
z (2s) = ip . 



^j 



Kj 



,-A.f 



Kj 



-Zj(s)=ipj, j = 1,2, 



□ 



4.2 Construction of symmetric periodic solutions 

Lemma T4 . 2 1 allows one to explicitly find all (s, s)-periodic solutions according to the following algorithm. 

Step 1. For each s > 0, we find the (unique) ipj — ^j( s ) such that Vj(s) = —ipj f° r J = 1,2,..., assuming 
that T-L(v) = 1 on the interval [0, s). To do so, we solve the equation (cf. Remark |2.6[) 

which yields 

K 1 — p _A i s 

^ = ^) = -a7-iT^- (4 ' 7) 

We note that ipj(0) = and ipj(s) monotonically decreases and tends to —Kj/Xj ass4 +oo. 
Step 2. We find the (unique) ipo = ipo(s) such that ip = a. To do so, we solve the equation 

oo 

m ipo + J~] rrijipj = a, 
i=i 

which yields 

1 / 00 

m^jis) I . (4.8) 



1 / °° 
ip = ip (s) = — a - V 



Note that the function ip with the Fourier coefficients given by (|4.7I) and (|4.8I) belongs to ff 1 . This follows 
from (|2~T4| and (|2~T6|) . 
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Step 3. If the solution^ v(x,t) = v(x,t; s) of problem (|2.9p - (l2.11[) with the initial data -0 = ip(s) is such 
that H(v) does not switch for t < s and switches at the moment t = s, then, by Lemma 14. 2 [ there exists a 
2s-periodic solution z(x, t; s) (which coincides with v{x, t; s) for t < s). 

The switching condition is 



^mjVj(s) = fi, 

or, equivalently, 



00 J r 2 —As 

F(s) := m Q K Q s + 2 V mj -r^ • - = Z 3 ~ ( 4 - 9 ) 

' A, 1 + e a j s 

3=1 J 



To check that the switching does not occur before s, we note that, due to Remark l2.6l the mean temperature 
v{t; s) corresponding to the initial condition (|4.7|) . (|4.8j) is given by 

°° °° K 1 - e~ A ^ 

Therefore, the condition s) = /3 is equivalent to 

m fc i + 2) m j — 1 - — t— = fi - a. 

Taking into account equality (|4.9p . we see that the condition v(t; s) = (3 is equivalent to the following: 

oo A s \ t 

H(t, s) := m k (t - s) + 2 V mj -^ ■ & ' ~_ F A ./ = 0. (4.10) 

\j l + e A i s 

3 = 1 J 

Moreover, the fulfillment of the inequality H (t, s) < for all t € (0, s) is necessary and sufficient for the absence 
of switching moments before the time moment s. 

Definition 4.2. We will say that F(s) and H (t, s) are the first and the second characteristic functions, 
while (|4.9|) and (|4.10[) are the /irsi and the second characteristic equations, respectively. 

The first and the second characteristic equations will play a fundamental role in the description of periodic 
solutions and their bifurcation sets (see Theorems 14.11 and 14.21 below) . 

The following lemmas describe some properties of the characteristic functions. 

Lemma 4.3. 1. F(s) is continuous for s > and analytic for s > 0, 

2. F(0) = 0, F(s) increases for all sufficiently large s > 0, and lim F(s) = +oo, 

s— >+oo 

3. for each fi — a > 0, the first characteristic equation (|4.9|) has finitely many roots, 

4- the positive zeroes of F(s) are isolated and may accumulate only at the origin. 

Proof. 1. The series in (|4.9p is absolutely and uniformly convergent for Res > due to the Cauchy- 
Bunyakovskii inequality and (I2.16[) . Therefore, F(s) is continuous for s > and analytic for s > 0. 
Assertion 2 is now straightforward. 

To prove assertion 3, we note that, for fi — a > 0, the (positive) roots of the first characteristic equation (14.91) 
cannot accumulate at the origin. This follows by the continuity and the relation F(0) = 0. The roots cannot 
accumulate at infinity either (due to the monotonicity for large s). Therefore, all the roots belong to a compact 
separated from the origin. Now the analyticity for s > implies assertion 3. 

Assertion 4 follows from the analyticity of F(s) for s > and from the monotonicity for large s. □ 

Similarly, one can prove the following lemma. 
Lemma 4.4. 1. H(t, s) is continuous for s > 0, < t < s, 

2. for each s > 0, H (t, s) is analytic in t for t > 0, 



1 Here and further, we sometimes write s after the semicolon to explicitly indicate that the function depends on the chosen first 
switching time s as on a parameter. 
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3. H(0, s) = -F(s) and H(s, s) = 0, 



4- if s > and F(s) > 0, then the second characteristic equation (|4.10l) has no more than finitely many 
roots in t for t G (0, s). 

Taking into account Lemmas 14.11 14.31 and 14.41 we formulate the above algorithm as the following theorem 
(also mind Remark |4.1[) . 

Theorem 4.1. 1. For a given ft — a > 0, there are no more than finitely many periodic solutions of prob- 
lem (|2~Tj) , (liOj) . which we denote z {1 \ . . . , z {N \ 

2. All the periodic solutions . . . , are symmetric. 

3. If si, . . . , sjv are half-periods of z^\ . . . , z^ N \ respectively, then Si, . . . , Sjv are the roots of the first char- 
acteristic equation (|4.9[) . 

4- Let sn+i, . . . , sjvi be positive roots of the first characteristic equation (I4.9P different from si, . . . , sn- TTien 

foj H(t, sj) < /or all t G (0, Sj-) ifj = l,...,N, 

(b) H(t; Sj) = /or some i e (0, sj) if j = N + 1, . . . , N%. 

In particular, Theorem l4.1l implies that a positive root Sj of the first characteristic equation (|4.9|) "generates" 
a 2s :) -periodic solution if and only if H (t, Sj) < for all t € (0, Sj). 

Now we will keep the domain Q and the functions m{x) and if (a:) fixed, while allow the thresholds a and ft 
vary. We will classify the existence of all periodic (i.e., (s, s)-periodic) solutions with respect to the parameter 
s and with respect to the parameter ft — a. By the existence of a periodic solution for a given s > we mean 
that there exist numbers a < ft (depending on s) such that problem (|2.1[) , (|2.3p with these a and /3 admits an 
(s, s)- or, equivalently, a 2s-periodic solution. 

First, we show that one can divide the positive s-semiaxis into intervals (whose union is denoted by L) in 
the following way. For every interval L C L, either there are no 2s-periodic solutions for all s € L' or there is 
exactly one 2s-periodic solution for every s G L', which smoothly depends on s m L' . The complement S of 
the union L of all those intervals will consist of points of possible bifurcation with respect to s (half-period). 
It will be a compact set. Typically, S will consist of finitely many points (see Examples 14. II) . 

The compact set £ = F(S) will consist of points of possible bifurcation with respect to the parameter ft — a. 
This set divides the positive (ft — a)-semiaxis into open intervals (whose union is denoted by A). For ft — a in 
an interval A' C A, the number of periodic solutions remains constant and they smoothly depend on ft — a G A' 
(see Example 14. lj) . 

First, we introduce the set 

So = {s > : F(s) = 0}. 

Due to Lemma 14.31 the set So consists of no more than countably many points, which may accumulate only at 
the origin. 

To introduce the next set, we denote for s > 

T(s)={te(0,s):H(t,s)=0}. (4.11) 

By Lemma 14.41 r(s) consists of finitely many roots of the equation H(t,s) = on the interval t G (0, s), 
provided that F(s) > 0. 
Consider the set 

Sx = {s > : F(s) > 0, t(s) = and H t (t,s)\ t=s = 0}. 

Thus, Si consists of those s for which the corresponding trajectory v(x,t; s) intersects the hyperplane tp — ft 
for the first time at the moment s and touches it nontransversally at this moment. Note that any number 
s G Si generates a 2s-periodic solution. 
Consider the set 

S 2 = {s > : F(s) > 0,t(s) ^ 0, and H t (t,s)\ t=t > = Vt' G r(s)}. 

Thus, S2 consists of those s for which the corresponding trajectory v(x,t; s) intersects the hyperplane cp = ft 
for the first time before the moment s and touches it nontransversally at each of the intersection moments 
(before s). None of the numbers s G S% generate a 2s-periodic solution. 
We also introduce the set 

S 3 = {s > : F(s) > and F'(s) = 0}. 
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We note that the set S3 consists of no more than countably many isolated points which may accumulate only 
at the origin. This follows from the analyticity of F(s) for s > and from the monotonicity for large s. 
Now we set 

L = (0, 00) \ S0US1 US 2 - 

and 

£ = F(Si U S 2 U S 3 ), A = (0,oo)\£. 

We note that the above sets Si, L and S, A do not depend on s or f3 — a. They only depend on rrij, Kj, and 
Xj. We also note that the sets Sq, . . .S3 and S are bounded. Indeed, Sq and S3 are bounded because F(s) 
monotonically increases for sufficiently large s. Furthermore, it is proved in |14) that, for sufficiently large 

dv(t; s) 

f3 — a (hence for sufficiently large s), the first switching moment for v(x, t; s) is equal to s and 



<it 

Therefore, S2 and S3 are also bounded. The boundedness of So, ■ ■ ■ , S3 implies the boundedness of £ 



> 0. 

t=s 



Theorem 4.2. 1. Let L 1 be an open interval in L. Then either there are no 2s-periodic solutions for all 
s G L' or, for any s G L 1 , there is a unique 2s-periodic solution z(x, t; s) of problem (12.1[) . (12.3[) . Moreover, 
the initial value z(x,0;s) smoothly depends on s G L' [in the H 1 -topology) . 



2. Let A' be an open interval in A. Then the number of periodic solutions of problem (|2.ip , (|2 . 3[) remains 
constant for all f3~a G A'. The initial values of those solutions and the first switching times continuously 
depend on (3 — a G A' (in the H 1 -topology) . 

Proof. 1. Let L' be an open interval in L. For any s G L', we denote by v(t;s) the mean temperature 
corresponding to the initial condition (|4.7[) , (|4.8p . We recall that 

v(t;s)=(3 

if and only if 

H(t,s) = 0. 

Fix an arbitrary s' G L' . Then s' £ S , i.e., F(s') ^ 0. If F(s') < 0, then F(s) < for all s G L' (otherwise, 
F(s) = for some s G L' , but then s G So, which is impossible). In this case, every s G L' does not generate a 
periodic solution. 

Assume that F(s) > 0. 

Consider the sets r(s) given by (|4. 1 1|) for s G L' . We claim that if t(s') = 0, then t(s) = in a sufficiently 
small neighborhood of s'; if t(s') 7^ 0, then r(s) 7^ in a sufficiently small neighborhood of s' . Indeed: 

la. Let t(s') = 0. Suppose that there is a sequence Si converging to s and a sequence ti G (0, Si) such that 
H(ti, Si) = 0. Taking a subsequence if needed, we can assume that U — > t' G (0, s']. Thus, by continuity 
of H(t, s), we have 

H(t',s')={) (4.12) 

Since t(s') — and s' ^ Si, we have H t (t, s')\ t=s > 7^ 0. Therefore, by the implicit function theorem and 
by the identity H(s,s) = 0, it follows that, in a neighborhood of the point (s',s'), the only root (in t) 
of the equation H(t,s) = is t = s. Hence, all U lie outside a fixed neighborhood of s', which means 
that t 1 < s'. Together with (|4.12p . this yields r(s') 7^ 0. This contradiction proves that r(s) — in a 
sufficiently small neighborhood of s'. 

lb. Now let r(s') ^ 0. Since s' S 2 , there is t' < s' such that H(t',s') = and H t (t, s')\ t = t > ^ 0. By the 
implicit function theorem the equation H(t, s) — admits a solution t = t(s) in a neighborhood of s' 
such that t' = t(s'). By regularity, t(s) < s if the neighborhood is small enough. Therefore, t(s) 7^ in 
a sufficiently small neighborhood of s' . 

To complete the proof of assertion 1, we choose an arbitrary compact interval in L' , cover each point of it 
by the above neighborhood and take a finite subcovering. 

The smooth dependence of the initial value of the periodic solution on s G L' follows from the explicit 
formulas (|4T71) and (|4T5|) . 

2. Let A' be an open interval in A. 

Fix an arbitrary b' G A'. Since b > 0, Lemma [4.31 implies that the first characteristic equation F(s) = b 1 
has finitely many (say, N{) positive roots s' 1; . . . , s' Ni . Since b' £ S, it follows that s[- £ S 3 , i.e., F'(s'.) 7^ 0. 
Therefore, for b in a neighborhood of b', there exist exactly Ni positive roots Si = Si(b), . . . ,sn 1 — SjVi(&) of 
the first characteristic equation F(s) — b, which smoothly depend on b. 
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Further, we assume that there are N (N < N\) numbers s[, . . . , s' N for which the minimal root of the 
equation H(t, s'j) — on the interval (0, s'j) is equal to s'j. As before, this means that s'j generate 2s^-periodic 
solutions for j = 1, . . . , N and do not generate periodic solutions for j = N + 1, . . . , iVi (cf. Theorem 14. ip . 

Since b' > and b' £ £, it follows that s'j $ S USi US 2 (j = l,...,Ni). Therefore, similarly to part 1 of the 
proof, for all b in a neighborhood of b' , the numbers sj — Sj(b) generate 2sj -periodic solutions for j — 1, . . . , N 
and do not generate periodic solutions for j = N + 1, . . . , N±. 

To complete the proof of assertion 2, we choose an arbitrary compact interval in A', cover each point b' of 
it by the above neighborhood and take a finite subcovering. □ 



Remark 4.2. Theorem 14.21 indicates the ways a new periodic solution may appear or an existing periodic 
solution may disappear, i.e., bifurcation occurs. 

When varying the parameter s, bifurcation may occur only if s G So U Si U S2. 

1. The condition s € Sq implies that a and (3 coalesce. 

2. The condition s G Si corresponds to the tangential approach of the trajectory v(x, t; s) to the hyperplane 
tp — P — a + F(s). At the point s, the periodic solution exists. In the literature on switching (or hybrid) 
systems, such a bifurcation is usually called "grazing bifurcation" . The corresponding Poincare map will 
be discontinuous at this point. 

3. The condition s £ S2 also corresponds to the tangential approach of the trajectory v(x,t;s) to the 
hyperplane ip = (3 = a + F(s). However, at the point s, the periodic solution does not exists. The 
switching occurs before the trajectory comes in the "symmetric" position. This bifurcation can also be 
called "grazing bifurcation" . 

When varying the parameter j5 — a > 0, bifurcation may occur if a point s G F~ 1 (/3 — a) belongs to Si, S2, 
or S 3 . 

Grazing bifurcation occurs on Si and S2 as described above. 

If s G S3 \ (Si U S2), then a new root of the first characteristic equation (|4.9|) may appear and then split 
into two roots (or two existing roots may merge into one and then disappear) as j3 — a crosses the value F(s). 
If the first switching moment for v(x,t;s) is equal to s (i.e., H(t,s) < for t < s or, equivalently, r(s) = 0), 
then a new periodic solution will appear and then split into two (or the two existing periodic solutions will 
merge into one and then disappear). This corresponds to a fold bifurcation. 

On the other hand, if the first switching moment for v{x, t] s) is less than s (i.e., H(t, s) = for some t < s 
or, equivalently, r(s) 7^ 0), then no bifurcation happens. 



We consider an example illustrating Theorems 14.11 and 14.21 

Example 4.1. Let Q be a one-dimensional domain, e.g., Q = (0,7r), cf. [8 HTT1|2"3"] . Let the boundary condi- 
tion (|2.3p be given by 

v x (0,t) = 0, v x (n,t) =H(v)(t). 

From the physical point of view, these boundary conditions model a thermocontrol process in a rod with 
heat-insulation on one end and a heating (cooling) element on the other. 
It is easy to find that 

A o =0, eo = ^, K = e (7T) = 

Xj=j 2 , ej(x) = J^cosjx, Kj = ej(Tr) = (-IV^, .7 = 1,2, .... 

Let mo = 2, mi = 7712 = 4, and m 3 = = • • • = 0. Then the bifurcation diagram is depicted in Fig. 14.11 
Let mo = 3.2, mi = m,2 = 4, and m.3 = 777,4 = • • • = 0. Then the bifurcation diagram is depicted in Fig. 14.21 
"Evolution" of periodic solutions with respect to the parameter j3 — a is visualized in Fig. 



In |14) . it was shown that there exists a unique periodic solution if f3 — a is large enough. Moreover, it is 
stable and is a global attractor. To conclude this section, we prove that a periodic solution can also exist for 
arbitrarily small /? — a. Further, we will show that such a solution need not be stable. 

Assume that the following condition holds. 

Condition 4.1. The functions m 6 H 1 and K € H 1 ' 2 satisfy 

00 „ 
M := V rrijKj = / m{x)K{x) dT > 0. 

JdQ 



3=0 
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trajectory corresponding trajectory corresponding 



to point A 



to point B 



P-a 



one periodic solution 



/3 - a « 0.23 
P - a ss 0.04 



two periodic solutions 
one periodic solution r 




'ghost" 
trajectory 



.4 




0.26 



4.10 



Figure 4.1: Bifurcation diagram for mo = 2, mi = r«2 = 4, and m3 = 7714 = • • • = 0. For any s > 0, there exists 
a unique 2s-periodic solution if the graph of F is bold at the point s and there are no 2s-periodic solutions 
otherwise. For any j3 — a > 0, there exist one or two periodic solutions depending on whether the horizontal 
line levelled at (3 — a intersects the bold part of the graph of F at one or two points, respectively. Point 
A (s « 0.26, /3 — a w 0.23) on the graph corresponds to s G Sx. There exists a corresponding 2s-periodic 
solution, whose trajectory is tangent to the hyperplane (p = /5 at the moment s (see the left inset). Point B 
(s w 4.10, (3 — a s» 0.04) on the graph corresponds to s £ S 2 ; there does not exist a 2s-periodic solution for this 
s. However, if one did not switch when v(t; s) tangentially intersected the hyperplane <p = j3 at the moment s, 
but switched only when v(t; s) intersected the hyperplane (p = (3 for the second time (at some moment Si > s), 
then the resulting trajectory would be 2si periodic. Such a trajectory is referred to as a "ghost" trajectory 
(see the right inset). 



The convergence of the sum follows from Remark l2.5l The equality follows from the definition of mj and Kj . 
The essential requirement of Condition 14. II is the positivity of the sum, or, equivalently, of the integral. From 
the physical viewpoint, this condition implies the presence of thermal sensors on a part of the boundary where 
the heating elements are. 

Theorem 4.3. Let Condition \4- l\ hold. Then there exist numbers lj > and a > such that, for any ft— a < uj, 
there exists a 2s-periodic solution z(x,t) = z(x,t;s) of problem (12.1[) . (|2.3[) such that s < a. On the interval 
(0,a;], the function s — s(/3 — a) is strictly monotonically increasing and s ~ > as /3 — a — > 0. 

00 

Proof. 1. By Condition 14.11 F'(0) = mjKj > 0. Therefore, for sufficiently small (3 — a > 0, the equation 

3=0 

F(s) = /3 — a has a unique solution s > in a small right-hand side neighborhood (0, a] of the origin. Clearly, 
the function s — s(/3 — a) possesses the properties from the theorem. 

Consider the solution v(x,t) — v(x,t;s) of problem (|2.1[) - (I2.3[) with the initial data ip — ip(s) defined in 
Steps 1-3 above. 

To complete the proof, it remains to show that v(t) = v(t; s) < (3 for t < s and apply Theorem 14. II 
2. Using representation (|2.19l) . Remark 12.61 and formulas (|4.7p . we have for t < s 

dv(t s) ^2, 2e" A J* ^2. / 2e~ X J t \ 

= m Ko + £ m j K j i + e _ XjS = M + ^ mjKj I — s -1 . (4.13) 

i=i j=l 

Using Remark [531 one can easily check that the absolute value of the series on the right-hand side is less than 
M/2 for sufficiently small s and t < s. Therefore, v(t; s) is monotonically increasing until the first switching 
moment. Thus, the first switching occurs for t = s. □ 



We stress that Theorem 14.31 ensures the uniqueness of a periodic solution with a small first switching time s 
(hence small f3 — a). However, the theorem does not forbid the existence of other periodic solutions with large 
period and large /3 — a. 
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Figure 4.2: Bifurcation diagram for mo = 3.2, mi = m-2 = 4, and 7713 = 7714 = • • • = 0. For any s > 0, 
there exists a unique 2s-periodic solution if the graph of F is bold at the point s and there are no 2s-periodic 
solutions otherwise. For any /3 — a > 0, there exist one, two, or three periodic solutions depending on whether 
the horizontal line levelled at j3 — a intersects the bold part of the graph of F at one, two, or three points, 
respectively. Points A (s ps 0.75, (3 — a ps 0.51) and B (s ps 1.74, /? — a ps 0.26) on the graph correspond to 
s £ S2. In each of these points, a 2s-periodic solution does not exist. However, if one did not switch when v(t; s) 
tangentially intersected the hyperplane <p = /3 at the moment s, but switched only when v(t; s) intersected the 
hyperplane ip — (3 for the second time (at some moment s± > s), then the resulting trajectory would be 
2s\ periodic. Such a trajectory is referred to as a "ghost" trajectory (see the insets). Point C (s ps 0.55, 
— a ps 0.56) corresponds to the fold bifurcation, where two periodic solutions merge into one and disappear 
as (3 — a increases and crosses the critical value ps 0.56. 

Remark 4.3. It is an open question whether one can choose the functions m{x) and K{x) and the parameters 
a and /3 in such a way that problem (|2.ip , (|2.3p has no periodic solutions. 

4.3 Stability of periodic solutions 

In this section, we will show that the thermocontrol problem with hysteresis may admit unstable periodic 
solutions. 

For simplicity, we assume that only finitely many Fourier coefficients rrij do not vanish (but see Remark 14. 7|) . 
Condition 4.2. There is N > 1 such that 

JJ = { rn , mi, . . .,TOjv}. 

Clearly, modifications needed if J consists of other Fourier coefficients rrij are trivial. 

Remark 4.4. The fulfilment of Condition 14.21 implies that m G H 1 . Moreover, the sum in Condition 14.11 
becomes finite: 

N 

^rrijKj > 0. 

Remark 4.5. If N = 0, i.e., J = {mo}, then it is easy to see that the (one-dimensional) guiding system (|2.2ip 
has a unique periodic solution for any a and j3 and this solution is uniformly exponentially stable. By Theo- 
rems 13.11 and 13.41 the same is true for the original problem (|2.ip , (|2.3I) . 
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Figure 4.3: Visualization of "evolution" of periodic solutions with respect to the parameter ft — a for mo = 3.2, 
m\ = mi = 4, and = = • ■ ■ = 0. For each /3 — a, the horizontal plane represents the phase space H 1 
with periodic solutions. Points A and B correspond to apparition (or termination) of periodic solutions, while 
point C corresponds to the fold bifurcation (cf. Fig. 14. 2} . 



Assume that Condition 14.21 holds. Let z(x,t) be a 2s-periodic solution of problem (|2.1[) . (|2.3[) . Denote by 
z(t) = (zo(t),z(t)) the corresponding 2s-periodic solution of the guiding system (|2.21l) . Let us study the map 
n a and the Poincare map II (see Sec. [3]) of the guiding system (12.21[) in a neighborhood of z(0). 

First of all, we consider the projections of these operators onto the TV-dimensional space V (see (|2.23ll ). 

We consider the orthogonal projector 

E: V -> V 

given by Eip = ip, where 

<p = {cpj}f =0 , <P = {<Pj}f=i- 
We also introduce the lifting operator 



R a : V -> V 



given by 



/ i N 

V m ° mo ti 



N 



Thus, R a E(y>) = ip for ip £ V such that m j ( Pj — a i an d ER a (^) = ip for cp G V (see Fig. 14.4 

3=0 

V 



t 




y E > 






-P^ ► 


— a 


= P 



Figure 4.4: The projection operator E and the lifting operators R Q and R,g 
Denote by II Q : V — > V the "projection" of n a onto V given by 

n Q (vO = En Q R a (¥j). 
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Similarly, one can define the operators R^ and lip. 

The operators E, R Q , and are continuously (and even inhnitcly) diffcrcntiable. Therefore, the opera- 
tors H a and Hp are also continuously differentiable, provided so are H a and 11,3 - 

We introduce the operator II : V — > V by the formula 

n( v ) = EnR Q ( v ). 

The following property of II is straightforward (see Fig. I4.5|) : 

n = UaU a . 







Km 


lip 

n(v>). 







{-Ki/Xu-Kj/Xj) 

Figure 4.5: The operators n a and II = n^n a in the space V — Span(ei, eg, . . . , ejv) 

It is easy to see that the point tp = z(0) is a fixed point of the map II acting in the iV-dimensional space V. 
In the formulation of the following results, we will use the following functions: 



Qj = QA s ) = l + e -\ jS > Q = ( 3( s ) = m ° K ° + m 3 K jQj( s ) 



,;=i 



We note that, due to (|4.13[) . we have at the switching moment s 

dz(t) 



dt 



= Q(s). 



(4.14) 



(4.15) 



In particular, this implies that Q(s) > 0. 

Theorem 4.4. Let Condition\4.£\ hold, and let z(x,t) be a 2s-periodic solution of problem (|2.1[) . (12. 3[) . Assume 



that Q(s) > 0. Then z(x,t) is stable [uniformly exponentially stable) if and only if the fixed point z(0) of the 
map II is so. 



Proof. Due to (|4.15l) . we have 



dz{t) 



dt 



„ dz(t) 

> 0. By symmetry, — - — 
t=s dt 



t=2s 



< 0. Now it remains to apply the 



formula IT(£>) = (R Q 1TE)(£>) and Theorem GEU □ 
To study the stability of the point tf> = z(0), we consider the derivative of II at the point ip. 

Lemma 4.5. Let Condition \4-2\ hold. If Q(s) > 0, then the operator H a : V — > V is differentiable in a 
neighborhood of xp = z(0) and the derivative 



at the point xp — z(0) is given by 

N 



N 1 / N 

D^IL a (ip)ip = Xl 'ipjej{x) + jrpr J]™ fc (l 

3 = l ^ [S) \fe=l 



N 



J2K J Q j (s)e j (x), 

) 3=1 



(4.16) 



where e%(x), . . . , ejv(x) form the basis in V and Qj(s) and Q(s) are defined in (|4.14[) . 
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Proof. Since Q(s) > 0, it follows from (|4.15p that 



dz(t) 



have 



dt 



> 0. Therefore, applying Lemma 4.2 in [2], we 



3 = 1 



dz(t) 



dt 



-1 / N 



N 



y^mfc (l 



\k=i 



i=i 



Taking into account equalities (I4.7[) . (|4.14p . and (|4.15p . we obtain the desired representation (|4.16[) . 



□ 



Remark 4.6. Due to Lemma T4.51 the linear operator D^H a (il)) is represented in the basis ei(x), . . . ,e^(x) 
by the (N x 7V)-matrix A = A(s) of the form 



(4.17) 



where 





(I 


- Ei + Sxo-\ 


S\o- 2 


Sia 3 


Sia N \ 






S 2 cri 1 - 


- E 2 + S 2 o~2 


S2V3 




A = 




S 3 ai 


5*30-2 1 - 


E 3 + S 3 a 3 


S 3 0~n 




\ 






Sn& 3 


1 — En + Sno~n ) 


Ej = 






Sj = Sj(s) = 


Q(s) ' 


°3 = Oj'O) = rnjEj{s) 



(4.18) 



Note that A(0) is the identity matrix. 

The following lemma results from Lemma 14.51 and from the symmetry of the periodic solution z(x,t). 



Lemma 4.6. Let Condition \4-^\ hold, and let Q(s) > 0. Then the operator II : V — > V is differ entiable in a 
neighborhood of tp — z(0) and the derivative 

A/,II(i/0 : V -> V 

at the point tp = z(0) is given in the basis ei(x), . . . , epj(x) by the matrix A 2 , where A is defined in (|4.17l) . 
Denote the eigenvalues of the matrix A = A(s) by /ii = A*i(s), i = 1, . . . ,N. 

The main result of this section is the following theorem. In particular, we will use it to construct unstable 
periodic solutions. 

Theorem 4.5. Let Condition \4.2\ hold, and letz(x,t) be a Is-periodic solution of problem (|2.1[) . (I2.3[) . Assume 
that Q(s) > 0. Then the following assertions are true. 

1. All the eigenvalues /Ltj of the matrix A satisfy m 7^ 1. 

2. If |^,| < 1 for all i — 1,...N, then the Is-periodic solution z{x,t) of problem ()2.1j) . (|2.3j) is uniformly 
exponentially stable. 



3. If there is an eigenvalue /ifc such that > 1, then the ^Is-periodic solution z(x,t) of problem (|2.1[) . (|2.3p 
is unstable. 



Proof. Assertion 1 follows from Lemma T4 . 71 below . It is known [16] that, under assumptions of items 2 and 3, 
a fixed point is, respectively, stable or unstable. By Theorem 14.41 and Lemma 14.61 this fact implies assertions 
2 and 3. □ 



Corollary 4.1. Let Conditions \ and \4-2\ hold. Then, for all sufficiently small (3 — a > 0, there exists a 
2s-periodic solution z{x,t) of problem (|2.ip , (|2.3p and assertions 1-3 in Theorem \4-5\ are true. 



Proof. The existence of z{x, t) follows from Theorem l4.3l Moreover, we have shown in the proof of Theorem[ 
that Q(s) > for all sufficiently small (3 — a > 0, provided that Condition 14.11 holds. Thus, the hypothesis of 
Theorem 14.51 are true. Therefore, the conclusions are also true. □ 



Now we prove the following auxiliary result, which we have already used in the proof of Theorem 14.5 
Lemma 4.7. Let Condition\4.2\ hold. IfQ(s) =/= 0, then the eigenvalues fii of A satisfy 



N 



N 



where Q is defined in (|4. 14|) and Ei in (|4. 18|) . 
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Proof. Substituting aj — rrijEj, we have 



N 



N 



j](^-i)=ia-ii =n^-iBi 



where 





/ 5*17711 — 1 


Sim 2 


Si m 3 


SlTTlN 


\ 




S 2 mi 


S 2 m 2 - 1 


527713 


S 2 rriN 




B = 


S 3 mi 


S 3 m 2 


5 3 m 3 - 1 . 


S 3 m N 






\ Sjymi 


SNm 2 


5jv™3 


■ S N m N - 


y 



and | • | stands for the determinant of a matrix. 
Let us compute the determinant of B: 





5*1 


S\m 2 


Sim 3 






s 2 


S 2 m 2 - 1 


S 2 m 3 


S 2 m N 


B| = mi 


s 3 


5 3 m 2 


S 3 m 3 - 1 . 


S 3 m N 






SNm 2 


S N m 3 





S 2 m 2 - 1 
S 3 m 2 

SNm 2 



S 2 m 3 
S 3 m 3 - 1 

S N m 3 



S 2 m.N 
S 3 m N 

Snitin — 1 



To find the determinant of the first matrix, we multiply its first column by rrij and subtract it from the jth 
column for all j = 2, . . . , N. As a result, we have 

S 2 m 2 - 1 S 2 m 3 . . . S 2 m,N 
S 3 m 2 S 3 m 3 - 1 ... 5*3771 



|B| = (-l^-^imr - 



SMm 2 



S N m 3 



Similarly decomposing the second determinant, we obtain (after finitely many steps) 

N m K 



|B| = (-lJ^OSiroi + • • • + S N m N - 1) = (-1) J 



Q 



□ 



Remark 4.7. Let us discuss modifications needed in the case of infinite set J in Condition l4.2l The construction 
of the maps n a , 11^,11 is quite similar and the modifications are obvious. The conclusion of Theorem 14.41 with 
the modified map II remains true. 

Formula (|4.16j) for the Frechet derivative D^,H a (if)) remains the same but the sums become infinite. Their 
convergence follows from Remark 12.41 Formally, the linear operator D^,H a (if)) can be represented as the 
matrix A (see (|4.17p ). which now becomes infinite-dimensional. 

It is proved in [14] that the operators Lla,]!^,]! are compact. Therefore, the same is true for their Frechet 
derivatives. In particular, this means that the spectrum of D^Tl a (xp) consists of no more than countably many 
eigenvalues, which may accumulate only at the origin. Thus, assertions 2 and 3 in Theorem 14.51 remain true 
(possibly with N = oo in assertion 2). 



4.4 Corollaries 

In this subsection, we assume that Condition 14.11 holds and that (3 — a > and s > are sufficiently small. 
Then Q(s) > and a 2s-periodic solution z(x,t) exists. Using Theorem 14. 51 we provide some explicit conditions 
of its stability or instability. Moreover, we will show that a periodic solution may have a saddle structure. 
The case N = is trivial (see Remark 14.51) . so we begin with the case N = 1. 



Corollary 4.2. Let Condition \4-2\ hold with N = 1. Then, for any (3 — a > 0, there exists a unique periodic 
solution z(x,t) of problem (12. lj) . (|2.3|) . The solution z{x,t) is uniformly exponentially stable. 



Proof. 1. By using the explicit formulas (Rermark 12. 6[) for the trajectories, we see that, for any trajectory 
v(x,t), the function v(t) either increases for all t > or first decreases and than increases. In particular, this 
implies that dv/dt > at the first switching moment. 

2. One can directly verify that the first characteristic function 

- Ai s 

F(s) := m Kos 



2mi — 
Ai 



1 



3 — XlS 



satisfies one of the two conditions: 



2G 



(a) F(s) > and increases for all s > 0, or 

(b) there is s* > such that F(s) < for < s < s* and F(s) > and increases for all s > s* . 

In both cases, the equation F(s) = f3 — a has exactly one positive root s\. 

3. Due to the observation in part 1 of the proof, the second characteristic function 

K\ e~ XlS - e~ Alt 

H(t, s) := m k (t - s) + 2m 1 — ■ — — 3— — = 

Xi 1 + e AlS 

satisfies the inequality H(t,S\) < for all t < si. Therefore, by Theorem 14.11 there is a unique 2si periodic 
solution of problem p7T]> . (f23|) . 

3. To prove its stability, we note that the matrix A consists of one element \x\. It satisfies (due to Lemma l4~7l 
or by direct computation) 

where Q > due to (|4.15|) and the observation in part 1 of the proof. If we show that /ii € (—1, 1), then the 
stability result will follow from Theorem 14.51 

Clearly, \i\ ^ \ for s\ > 0. One can also show that /ii ^ —1 for si > 0. To do so, one can check for 
example that the equation fii = — 1 uniquely determines m\K\ as a function of the other parameters. Then 
substituting it into the formula for F{s\) yields the contradiction F{s\) < 0. 

Since //i ^ ±1, fix G (—1,1) for sufficiently large si, and \i\ continuously depends on s%, it follows that 
jttx G (—1, 1) for any s±. □ 

Now we consider the case N — 2. 
Corollary 4.3. Lei Condition \4-2\ hold with N — 2, and let 

M = m K + m x K x + m 2 K 2 > 0. (4.19) 



Then, for all sufficiently small f3 — a > 0, there exists a 2s-periodic solution z{x,t) of problem (12. ip . (|2.3p 
uniquely determined by Theorem \4-3\ If 

(M - mi#i)Ai + (M - m 2 K 2 )X 2 < 0, (4.20) 
then I/X2I > 1 and z{x,t) is unstable for all sufficiently small f3 — a > 0. J/ 

(M - m 1 K 1 )\ 1 + (M - m 2 K 2 )X 2 > 0, (4.21) 
then \fii\, \^2\ < 1 a^rf z(x,t) is exponentially stable for all sufficiently small (3 ~ a > 0. 

Proof. 1. The matrix A is a (2 x 2)-matrix. Therefore, it has two eigenvalues /ii and fi 2 , which are either both 

real or complex conjugate. Denote S = 6i_ 2 = [i\. 2 — 1. Clearly, 8\^ 2 are the eigenvalues of A — I; hence, they 
are the roots of the quadratic equation 

(5 2 -tr(A-I)(5 + |A-I| =0. (4.22) 
Let us compute tr (A - I) and |A - I|. Due to (f47TT|) and (|4TT8]l . 

ll(A _ I) . El (_ 1 + !^) +E2 (_ 1 + ^ 

On the other hand, formulas (|4.18l) and (|4.14[) imply that Ej = Xjs + 0(s 2 ), and Qj(s) = 1 + 0(a), and 
Q{s) = M + 0(a). Therefore, 

tr (A-I) = a(\ 1+ 0(a)) (-1 + ^ + 0(s)\ +a(X 2 +0(a)) (-1 + ^ + 0(sfj = -(Ls+0(a 2 )), (4.23) 
where 

L = M _1 ((M - m 1 K 1 )X 1 + (M - m 2 K 2 )X 2 ). 

Further, by Lemma 14.71 



EiE2^- = s 2 (X 1 + 0(s))(X 2 + 0(a)) l^-j^ + 0(a) ) = J's' + 0(a% 14.24 ) 
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where 

J' = AiA 



2 * , m K 



M 

It follows from (|423|) and (l4~24l) that Eq. (|422|) is equivalent to the following: 

S 2 + (Ls + 0(s 2 ))S + J 2 s 2 + 0{s 3 ) = 0. 



Thus, 



Ls Vg - 4 J 2 + Q(s) 



2. If inequality (14.20j) holds, then L < and Re/ii, 2 > 1 for all small s > 0. 

Assume that inequality (|4.21[) holds, i.e., L > 0. If L 2 — 4J 2 + O(s) > 0, then the eigenvalues /xi^ are real 
and belong to the interval (0, 1). If L 2 — 4J 2 + O(s) < 0, then /ix,2 are complex conjugate and 

(Re iii) 2 + (Im/ii) 2 = 1 - Ls + (9(s 2 ) < 1, 

i.e., |/4i 2 1 < !• □ 



Example 4.2. Consider the problem described in Example 14.11 

Let mo > 0, mi = m-2 > 0, and m.3 — tua = ■ ■ ■ =0. Then condition (|4.19j) holds. Therefore, condi- 
tion (|4.20p . which implies the instability of the periodic solution for small s, takes the form 

mo < /^ A 2 - Ai _ 3\/2 
mi A 2 + Ai 5 ' 

while condition (|4.21l) . which implies the uniform exponential stability of the periodic solution for small s, takes 
the form _ 

wq > - Ai _ 3v/2 

mi A2 + Ai 5 

Finally, we show that periodic solutions can be unstable for N > 3. Moreover, if N is odd, they may have 
a saddle structure. 

Corollary 4.4. Let Condition \4-S\ hold with N > 3, and Ze£ 

N 

M = J2 m 3 K 3 > °- ( 4 - 25 ) 
J=0 



Then, for all sufficiently small j5 — a > 0, £/iere errasfc a 2s-periodic solution z(x,t) of problem (|2.1[) . (|2.3I) 
uniquely determined by Theorem ^. 3\ If 

N 

J2(M - mjK^Xj < 0, (4.26) 
j=i 

£/ien z(a;,t) is unstable. 

If we additionally assume that N is odd, then there is an eigenvalue of Z?^,II(z(0)) with real part greater 
than 1 and a real eigenvalue in the interval (0, 1). 

Proof. 1. The matrix A is an (N x A)-matrix. Due to 633), (|4TT8]) . and (t4T26|) . 

JV AT 

^ /*,- = tr A = AT + ^(5 3 m, - l)E, 
3=1 3=1 

N 

= N - M~ l 5^(M - mjK^XjS + 0(s 2 ) > N 

for sufficiently small s > 0. Therefore, the real part of at least one eigenvalue is greater than 1. By Theorem l4.5[ 
this implies the instability of z(x,t). 

2. Now we additionally assume that N is odd. By Lemma 14771 

JV 

no* - 1) < 0. 

3=1 
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Since iV is odd, the set of eigenvalues of A consists of an odd number of real eigenvalues . . . , fj,L (1 < L < N) 
and (N — L)/2 pairs of complex conjugate eigenvalues. Therefore, 

L 

JJO*i-i)<o. 

3=1 

Hence, there is at least one eigenvalue, e.g., fj,i, which is real and is less than 1. Taking into account that 
jLtj(O) = 1 and fij(s) continuously depend on s, we see that fi\ £ (0, 1). Applying Lemma 14.61 we complete the 
proof. □ 
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